ACEEE Int. J. on Information Technology, Vol. 02, No. 02, April 2012 



Automatic Relative Radiometric Normalization for 
Change Detection of Satellite Imagery 

Munmun Baisantry* 1 , Dr.D.S.Negi* 2 , O.P.Manocha* 3 
"Image Analysis Centre, Defence Electronics Applications Lab 
Defence Research and Development Organisation, Raipur Road, Dehradun-248001 
munmunbaisantry@gmail.com 



Abstract — Several relative radiometric normalization (RRN) 
techniques have been proposed till date most of which involve 
selection of pseudo invariant features whose reflectance are 
nearly invariant from image to image and are independent of 
seasonal cycles. Extraction of such points is quiet tedious and 
human operator has to provide mutual correspondence by 
choosing easily recognizable and time invariant points. In 
this paper, we intend to propose a new automatic radiometric 
normalization technique to select PIFs in panchromatic 
images known as Bin-Division Method. For multispectral 
images, MAD (Multivariate Alteration Detection) has been 
employed for selecting PIFs based on the assumption that 
MAD components are invariant to affine transformation. This, 
followed by robust linear regression constitutes the whole 
automatic radiometric normalization procedure. 

Index Terms — Pseudo-invariant features, Bin-Division 
Method, Canonical Correlation Analysis, MAD 
Transformation, Linear Regression, RMS Error criterion. 

I. Introduction 

Radiometric differences introduced in multi-date images 
make it extremely difficult to assess actual changes as exuded 
by changes in spectral reflectance, which in turn, makes 
radiometric normalization a vital step in satellite image 
processing process like change detection and analysis. 
Radiometric normalization can be categorized as absolute and 
relative. Absolute normalization, though sophisticated and 
accurate, requires sensor calibration coefficients, an 
atmospheric correction algorithm and related input data like 
atmospheric measurements at the time of image acquisition 
which are very difficult to avail. This shifts our focus to RRN 
which involves normalizing the brightness values of the 
image, band-by-band to the reference image of a different 
acquisition date such that the two images seem to be acquired 
with the same sensor under similar atmospheric conditions. 
[1]. Non-linear correction, method like Histogram Matching 
is generally not preferred as they produce a lot of spectral 
variations which can be confused as actual changes [2]. 

Keeping this premise in mind, the normalization of subject 
image can be carried out as: 

Yf>(x.?)=o-g*Ts(x : y) (1) 
Where F' (x, y) the normalized gray value of subject image is, 
F (x, y) is the actual gray value and g,o are the normalization 
coefficients namely gains and offset values. Several relative 
radiometric normalization techniques like Pseudo-invariant 
features [3], Radiometri Control set [4], No Change set using 
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scatter gram [5] have been proposed till date. These methods 
involve radiometric stable features known as Pseudo- 
invariant points to approximate the DN values. These features 
have same reflectivity even during different seasonal cycles 
and their brightness values are linearly distributed [3] . The 
Dataset used. The image dataset used in this study is obtained 
from IRS- ID satellite. The dataset comprises of two 
geometrically registered panchromatic images of same area 
taken on two different dates. For normalization of multispectral 
images, two geometrically registered images of another area 
taken on two different dates have been used. Fig. 3 (a) & (b) 
as well as Fig.4 (a) & (b) show the dataset used. 

III. Methodology 

The whole process of relative radiometric normalization 
can be divided into three steps: (i) pre-processing (geometric 
registration); (ii) selection of pseudo-invariant features (iii) 
determination of normalization coefficients from selected PIFs. 

A Pre-processing: Geometric registration 

Panchromatic and multi-spectral subject images are 
geometrically registered with the respective reference images 
before automatic relative radiometric normalization process. 

B. Selection of Pseudo-invariant Features 

In this article it is assumed that linear effects are much 
dominant than nonlinear and PIFs are those features in an 
image which are statistically invariant to any seasonal cycles; 
their reflectivity is almost linearly distributed in nature, either 
increasing or decreasing. They can either be bright set like 
urban features, concrete, rock structures etc or dark set like 
deep water bodies. [3] Following sections provide the 
description of proposed different approaches for pseudo 
invariant features: 

1) Bin-Division Approach: Bin-Division is an effective 
method to separate PIFs features from a set of pixels in 
panchromatic images. It assumes the fact that radiance 
reaching at sensor can be related as linear function of 
reflectivity of target on ground. Using these assumptions, 
the pixels of subject and reference images are divided into 
three categories. This is based on the brightness value of a 
pixel of subject image which has either increased, decreased 
or has not changed with respect to the brightness value of 
the pixel at the same spatial position of the reference image. 
The corresponding pixels of both the images are placed in 
three bins. Let Fs(x, y) and Fr(x,y) are subject and reference 
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images respectively. The bin divisions for both images pixels 
are as follows: 



F=(x : y)&FR(x : ,y) 
F 2.(t : y)<F E (x : y) 
F s (x : y)F=F E (x : .y) 



(2) 
(3) 
(4) 



The pixels described by equation 2,3,4 are bin 1 , bin2 and 
bin3 respectively. The pixels which exist in binl and bin2 are 
combination of pixels which have experienced genuine land 
cover changes in their brightness values and radiometric 
changes. So binl is union of positively changed pixels and 
positive linear PIFs pixels and bin2 is union of negatively 
changed pixels and negative linear PIFs. The proposed Bin- 
Division methodology separate out PIFs from changed pixels 
in binl and bin2 assuming that the PIFs are linearly 
distributed. In Bin-Division method, our aim is to compute 
PIFs by minimizing the cost function that gives the linear 
similarity for corresponding pixels in bins. Let f (x,y,n) and 
f (x,yn) are pixels in bins and f (x,y,n) are normalized values, 

then cost function is defined as Cost (1,2,3 n) = Dist 

[f.(x,y,n) , f (x,y,n) , f (x,y,n) , W(o,g)] for n values of bin. The 

PIFs are the pairs for which Cost( 1,2,3 n)<a, a a user 

defined threshold that sets the level of normalization. The 
minimization of cost function depends upon 
Dist[f,f ,f ,W(o,g)] and weight W. The following section 
describes the determination of weight function algebraically. 
The objective of this method is to estimate the weights 
function W, for given f (x,y,n) and f (x,y,n), some knowledge 
of cost function. By assuming these quantities, it is possible 
to formulate a unified linear algebraic framework. The central 
of algebraic approach is the concept of finding a good 
estimate of weight function W which minimizes the predefined 
criterion of cost function defined by: f-f(W) =C. 

In the absence of any knowledge of C, our need is to find 
a criterion function to find W such that f (W) approximates f 
by assumption that norm of C is small as possible. In other 
words we wish to find out W such that criterion function 
IICIP=II f - f (W) IPis minimized. The minimizing criterion function 
isJ(W) =11 f- f (W) IP with respect W. Minimizing this equation 

we get £ ™=0 = -2/;</ r -/,(lV)} .Solving this 



equation, we get 



(5) 



The procedure is applied for each for 'Greater than' binl 
as well as 'less than' bin2 and pixels have residuals greater 
than a user-defined threshold were filtered out, leaving 
behind a final set of PIFs for both the bins. Closer the residual 
closer to the threshold value of 1, better would be PIFs. The 
plots shown in Fig. 1 (a) and (b) show the linear nature of the 
final PIFs selected from both the bins. 




-* * 4 * *- 

Fig.l.(a) A plot of PIFs selected for negative bin 




Fig.l.(b) A plot of PIFs selected by positive bin 

Normalization coefficients are calculated for both the bins 

from the resultant PIF features. From equation 5 [9], 
normalization coefficients for both the bins are as follows: 



o T = ■ 



o = 



(?) 
(8) 
(9) 



Where o + ,o ,g + ,g ,f + ,f ,f + ,f ,n ,n" are ordinates, gains, gray 
level values of subject, reference pixels and total number of 
pixels for respective bins. Using normalization parameters, 
normalized image is generated. 

Bin-Division method works well even when temporal 
variations are significant and relies on the pixels function, 
which excludes pixels marked with spectral changes. 

2) Multivariate Alteration Detection (MAD): In MAD 
automatic selection of PIFs is based upon the assumption 
that atmospheric and calibration differences are linearly related 
and MAD components are invariant to such affine 
transformations [7], [8]. This method is based on Canonical 
Correlation Analysis (CCA) [8]. 

Assuming [7] that we have two temporal images having N 
spectral channels, we can represent intensities of both images 
by random vectors F and G, where 

F=[F p F,F f F N fand G=[G,,G 2 ,G 3 G N f 
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Combining the intensities from all the bands linearly, we 
achieve 

U=aT=[ ai F 1+ a 2 F 2 + a 3 F 3 + a N F N ] (10) 



V=b T G= [b 1 G 1 +b 2 Gvfb 3 G 3 + b N G N ] 



(11) 



The differences of the linear combinations U and V are 
known as MAD variates. Maximizing the difference 
information requires maximizing variance which requires 
minimizing the correlation between U and V under the 
Var(U)=Var(V)=l with keeping the correlation as positive. 
Using CCA, searches for linear combinations U=a T F and 
V=b T G of the (ideally) Gaussian distributed variables 
[U,Vl€N(|i, X . ) with maximum correlation: 

(12] 



p=CoMu,v] =- 



This leads to the coupled generalized Eigen value problems 



^f3^S3^Sf a - P ^ff a 



(13) 

(14) 



MAD Transformation is defined as the subtraction of 
canonical variates (solutions of (13) and (14)) . The MAD 
variates are uncorrelated to each other, the last component 
being maximally correlated with minimum difference 
information and first component having the maximum 
difference information. 

Two images acquired under different atmospheric 
conditions with no genuine land cover differences will have 
a difference component with a nearly Gaussian distribution. 
Thus MAD components, in such cases have an uncorrelated 
Gaussian distribution. In case some land cover changes have 
occurred, the distribution will deviate more or less from its 
Gaussian nature. 

Using suitable decision threshold, the non-difference 
pixels can be separated from changed pixels. Furthermore, [8] 
the sum of squared standardized MAD variates is chi-square 
distributed, which helps in setting the threshold as: 



— v 



(15) 



where t=± 2 p=00l , with P as the probability of observing that 
value of t or lower. As long as the radiometric effects are 
linear, the pixels chosen can serve as PIFs. The plots shown 
in Fig. 2(a) and (b) depict nearly linear nature of PIFs for each 
band. The corresponding pixels values which are less than 
the threshold t, are PIFs for respective band. 
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Fig.2.(a) A plot of PIFs for Bandl 
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Fig.2.(b) A plot of PIFs for 



Band2 




Fig.2.(c) A plot of PIFs for 



■I* 

Band3 



rV. RESULTS AND DISCUSSION 

The quality of radiometric normalization in both the cases 
is statistically assessed using Root Mean Square Error 
(RMSE) which is defined as: 



RMSE k 



(16) 



Fig-3(a)-(c) and Fig-4(a)-(c), show panchromatic reference, 
subject, normalized subject images as well as MSS reference, 
subject, normalized subject images respectively. From Table- 
3 and 4, the RMS errors between the respective images( pan 
and mss) using proposed methods and other conventional 
methods can be seen. Numerically reduced RMSE after 
normalization in Table-3 One can observe that the RMS errors 
for proposed techniques are less as compared to other 
methods. 

Table I: Normalization coefficients for panchromatic image 







o(ordinite; 


Positive 


10.997 


:.cc: 


Xezative 


-11 512 


0.954 



Table II: Normalization coefficients for multispectral image 



BandXo. 


g(gaiiO 


ciordtuste; 


Band 1 


1.776 


-10S.329 


Band 2 


:.:ss 


25.605 


Band 3 


1251 


10114 



Table III: RSME using Bin-Division, Histogram Matching and Simple 
Regression 





Root Mean Square Error 


Bsfo-r; XonnalizatiMi 


1436 


ycrni'lizaicn (Bin-Division) 


7.72 


Ncmisliziticn (Histogram Matching) 


16.2923 


Xormalization (Simple Regression) 


9.6S 
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Table IV: RSME using MAD and Histogram Matching 



S.No 


Before 

Normalization 


Afta 
Normalization 
Histogram 
Matching 


After 

Normaliiation 
(Using MAD) 


Bandl 


13 04 


12.2 


1190 


Band2 


31.99 


21.79 




Band 3 


2S.S0 


1963 


176 


Ave. 


2466 


17.S7 


1632 





Eg 3(a): Reference MSS Image Fig 3(b): Subject MSS Image 




Fig 3(c): Normalized MSS Image 




Fig 4(a): Reference PAN Image 



Fig 4(b): Subject PAN Image 



Fig 4(c): Normalized PAN Image 



V. CONCLUSION 

The proposed techniques for radiometric normalization 
are based on the assumption that atmospheric and calibration 
differences are linearly related and pseudo-invariant features 
are invariant with time or linearly increasing or decreasing 

The techniques are computationally efficient and faster 
and are capable of radiometrically normalizing panchromatic 
as well as multi spectral images automatically with least human 
intervention, except for selection of threshold and it also 
gives distributed sufficient number of PIFs over the 
radiometric range of image. Statistical assessment tools like 
jRJVISE as well as visual interpretation prove that both the 
techniques can be developed as powerful framework in 
removing radiometric differences between two images, making 
them more visually similar. Comparison with RMSE for 
Histogram Matching and Simple Regression prove our results 
further. 
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